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Abstract Theoretical predictions for scattering processes with multi-particle final states at 
next-to-leading order (NLO) in perturbative QCD are essential to fully exploit the physics 
potential of present and future high-energy colliders. The status of NLO QCD calculations and 
• (-( ■. tools is reviewed. 

&1 1 Introduction 

The study of hard scattering processes at the Large Hadron Collider (LHC) pLj and a future 
TeV-scale linear collider is our primary means to probe and extend the Standard Model of 
particle physics. It is driven by the comparison of experimental measurements with theoretical 
predictions, which depends on our ability to compute collider cross sections in perturbative 
QCD with adequate accuracy [21 [3]. This can only be achieved by going beyond leading order 
(LO) in QCD. When using conventional measures, LO scale uncertainties are typically large 
compared to experimental uncertainties. Moreover, for theoretical reasons a reliable estimation 
of the scale uncertainty is not feasible at LO. Consequently, an assessment of different scale 
choices, which is particularly important for many-particle/jet processes, is not possible. Fur- 
thermore, the convergence of the perturbative series cannot be assessed at LO. When going 
beyond LO by including NLO corrections, the situation improves significantly!! At NLO, scale 
uncertainties can be assessed more reliably, and the residual uncertainties are often comparable 
to experimental uncertainties!! NLO calculations thus deliver accurate predictions not only 
for the overall normalisation, but also for kinematic distributions including peripheral phase 
space regions. This is in part due to the fact that new subprocesses often become active at 
NLO, which modify the normalisation and kinematic distributions. Our ability to determine 
the uncertainty of parton distribution functions (PDF) and to model the structure of jets is 
also greatly enhanced at NLO. 



^Presented at Linear Collider 2011: Understanding QCD at Linear Colliders in searching for old and new 
physics, 12-16 September 2011, ECT*, Trento, Italy 

^For processes with vastly differing scales, the resummation of large logarithms of ratios of scales may also 
be necessary. 

^ Notable exceptions are the hadroproduction of Higgs and Whh with (Tnlg/o'lo ~ 2. 



1 



LCll Proceedings 



Frascati Physics Series - Vol. LIII 



In Section [2], the state-of-the-art methods, implementations and tools for parton-level NLO 
calculations are briefly reviewed. In Section |3l the status of collider physics applications is 
described. The review ends with a summary]^ 



2 Methods, implementations and tools 

The structure and imphed modularity of NLO calculations is illustrated in Eqs. (l)-(3): 



CNLO — CBorn + 0" corr (1) 
(XBorn = j d(j)n ^ |AoP (2) 

(^corr = / d(f)n ^ |^ / d(j)jVj + ^LO-^nlo.V + -^LO^^NLO.V 

+ I g [|-Mnlo,rP - E^^] (3) 

The new components of the NLO correction acorr arejf] the virtual corrections (involving one- 
loop amplitudes), the real corrections (involving tree amplitudes) and the infrared subtraction 
termsj^ The resulting procedure for NLO calculations is given in Table [H The Binoth Les 
Houches Accord, a standard interface for combining the tree-level and loop-level contributions, 
has been defined in Ref. j6] and is implemented in many automated tools (see below). 

Until circa 2005, the limiting factor of NLO calculations was the computation of the virtual 
corrections, which typically applied Passarino-Veltman (PV) [7] or PV-inspired [8] tensor in- 
tegral reduction methods to evaluate the form factors of a Feynman-diagram-based amplitude 
representation. Several one-loop integral libraries are available as public codes: Loop Tools 
[TO], QCDLoop [n], Golem95 [12], OneLOop [I3] and PJFry [11]. The PV approach is 
general, but practical limitations arise due to the factorial growth of the number of Feynman 
graphs with N = n + 2, the strong growth of the number of reduction terms with and due to 
numerical instabilities for exceptional kinematic configurations, which are caused by vanishing 
Gram determinants. It has nevertheless been used successfully to create collections of NLO 
calculations based on analytic formulae and semi-automated methods, such as MCFM [T5| [T6]. 
MC@NLO [H] and VBFNLO [HI HH |20l |2l] Since 2004, tremendous improvements have been 
achieved for the calculation of multi-leg one-loop amplitudes due to the exploitation of on-shell 



^The important topics of next-to-next-to-leading order (NNLO) calculations and combining parton-level 
fixed-order calculations and parton-shower event generators are beyond the scope of this review. 
^The Born amplitude is assumed to be at tree level. 

^An alternative to the widely used subtraction formalism [4] is the phase space slicing method [5]. 
''The POWHEG BOX 22 library project [23l [24] was inspired by these collections. 
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1. Real correction: generate and evaluate 2 — )■ n + 1 tree-level amplitudes 

2. Subtract soft and coUinear singularities due to single unresolved real radiation 
to obtain finite result 

3. Integrate over (n + l)-particle phase space 

4. Virtual correction: generate and evaluate UV-renormalised 2 — )■ n one-loop 
amplitude after extraction of soft and collinear singularities to obtain finite 
result 

5. Confirm cancellation of soft /collinear singularities (absorb initial state collinear 
singularities into PDF) 

6. Integrate over n-particle phase space 

7. Combine 2 — >■ n -|- 1 and 2 — )■ n contributions 

8. Convolve with NLO PDF 

9. Repeat for all contributing subprocesses 



Table 1: Steps to calculate the NLO QCD corrections for a 2 — n process, n excludes elec- 
troweak decays. 



recursion relations and generalized-unitarity-cut constructibility as well as the possibility to 
even reconstruct the full rational terms |25l [26] . On-shell reduction related tools are Cut Tools 
[27], Rocket [28] and Samurai [29]. Further innovative, complementary methods are also being 
developed [30]. A comprehensive review of methods for multi-leg one-loop calculations can be 
found in Ref. [31] . 

Three widely-used algorithms for the generation of process-independent infrared subtraction 
terms are Catani-Seymour dipole subtraction [32], Frixione-Kunszt- Signer (FKS) subtraction 
[33] and antenna subtraction [M]]^ Several implementations for these standard schemes are 
available: Sherpa-Dipoles [36], MadDipole [37], HELAC-Dipoles [38j, MadFKS [39], TeVJet 
[10] and AutoDipole [5T]. 

The following programs aim to provide a comprehensive, automated solution for NLO cal- 
culations: aMC@NLO [211 EE 112], BlackHat/Sherpa [2S1 ESI US] , HELAC-NLO [lEl EZl EHl IB] , 
GoSam [15], FeynArts/FormCalc/LoopTools [101116] and MadGolem [17]. 

3 Collider physics applications 

Discussions at the Les Houches 2005 Physics at TeV Colliders Workshop resulted in a list of 
processes for which the knowledge of NLO corrections was considered of particular importance 

* Research on alternative subtraction schemes is also being carried out [35] , 
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for the LHC physics programme [IH]. This experimenter's NLO "wish hst" has guided theo- 
retical efforts and was subsequently revised and updated in 2007 |19] as well as 2009 [50j . The 
most recent version is displayed in Table [2j 

Due to the groundbreaking advances outlined in Section [2|, since 2009 the frontier for collider 
physics applications of NLO techniques has also advanced considerably. The following 2 — )■ 4 
processes - most are on the wish list - have now been calculated at NLO 

[21], ^ jets [621E31E61EI], ^ ^, 7* +3 jets [68] , ^ [591 EQl iB [69] , ^ ttjj 
[61 [70], PP ^ bbbb [71], PP ^ W+W-bb [72], ^ W^W^jj [21 [73], pp ^ W+W-jj [H] 
and most recently pp — )■ 4 jets [75]. Leptonic decays of weak bosons can be included trivially. 
At the same level of complexity, complete off-shell effects for pp — t- ti with dileptonic decay, i.e. 
pp — )■ e~^i'ebfi~Ufj,b, have been calculated at NLO QCD in Ref. [76j, which allowed to explicitly 
confirm the 0{asT/M) effect predicted by Ref. [77]. Advancing the frontier for linear collider 
physics, the process e~^e~ — )■ 5 jets has recently been calculated at NLO [78], which allowed 
to extract a competitive value of as{Mz) from 5-jet LEP data. Going beyond 4-particle final 
states in general requires the computation of 7-point one-loop amplitudes or higher. This is the 
current complexity frontier. At this level, NLO cross sections in leading-colour approximation 
have been calculated for 1^ + 4 jets by the BlackHat/Sherpa collaboration {pp — )■ 1^ + 4 jets [75] 
and pp ^ Z + 4: jets [80]) and for e~^e~ — >■ n jets up to n = 7 [81] The n = 7 case required 
the computation of a one-loop 8-point function. 



4 Summary 

NLO QCD predictions for multi-particle processes are essential to fully exploit the physics 
potential of the LHC and a future linear collider. In recent years, tremendous progress has 
been made in developing the calculational methods and tools that are required to compute 
NLO corrections for hard scattering processes with 6, 7 or more external particles. At this level 
a (semi-)manual approach is no longer feasible, and the transition from collections of codes for 
specific processes to automated code generation for any process up to a maximum complexity 
has now been achieved. Several such automated tools are available or will become public in the 
near future. The modularity of NLO calculations allows to interface many tool components on 
the basis of the Binoth Les Houches Accord. 

^pp is given as initial state, but pp is also implied. 

^"Recently, the full-colour virtual contribution to pp ^ W + 4 jets has been calculated [5^ . 
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Process {V £ {Z,W,-/}) 


Comments 


Calculations completed since Les Houches 2005 




1. -> VV+jet 

2. pp Higgs+2jets 

3. pp^VVV 

4. pp tt bb 

5. pp y+3jets 


H^VF+jet completed by Dittmaier/Kallweit/Uwer [511 [5^: 
Campbell/Ellis/Zanderighi [53]. 
ZZ+jet completed by 

Binoth/ Gleisberg/Karg/Kauer/Sanguinetti [54] 

NLO QCD to the gg channel 

completed by Campbell/Ellis/Zanderighi [16j: 

NLO QCD+EW to the VBF channel 

completed by Ciccolini/Denner/Dittmaier [55l [56] 

ZZZ completed by Lazopoulos/Melnikov/Petriello [57] 

and WWZ by Hankele/Zeppenfeld H9j 

(see also Binoth/Ossola/Papadopoulos/Pittau [58]) 

relevant for ttH computed by 

Bredenstein/Denner/Dittmaier/Pozzorini [59] [60] 

and Bevilacqua/Czakon/Papadopoulos/Pittau/Worek [61] 

calculated by the Blackhat/Sherpa [62] 

and Rocket ^63j collaborations 


Calculations remaining from Les Houches 2005 




6. pp ti+2jets 

7. pp -> VVbb, 

8. pp V^y+2jets 


relevant for ttH computed by 
Bevilacqua/Czakon/Papadopoulos/Worek [64] 
relevant for VBF H VV , ttH 
relevant for VBF H ^VV 
VBF contributions calculated by 
(Bozzi/) Jiiger / Oleari / Zeppenfeld [201 


NLO calculations added to list in 2007 




9. pp ^ bbbb 


qq channel calculated by Golem collaboration [65] 


NLO calculations added to list in 2009 




10. pp y+4jets 

11. pp Wbbj 

12. pp titi 


top pair production, various new physics signatures 
top, new physics signatures 
various new physics signatures 


Calculations beyond NLO added in 2007 




13. gg W*W* 0{a^al) 

14. NNLO pp ^ ti 

15. NNLO to VBF and Z/7+jet 


backgrounds to Higgs 

normalisation of a benchmark process 

Higgs couplings and SM benchmark 


Calculations including electroweak effects 




16. NNLO QCD+NLO EW for W/Z 


precision calculation of a SM benchmark 



Table 2: The experimenter's wish list for LHC processes in early 2010 (from [50j). 
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